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Abstract. In our previous work [20], [H] we introduced a correc- 
tion to the mean field approximation of interacting Bosons. This 
correction describes the evolution of pairs of particles that leave 
the condensate and subsequently evolve on a background formed 
by the condensate. In we carried out the analysis assuming 
that the interactions are independent of the number of particles 
N. Here we consider the case of stronger interactions. We offer 
a new transparent derivation for the evolution of pair excitations. 
Indeed, we obtain a pair of linear equations describing their evo- 
lution. Furthermore, we obtain apriory estimates independent of 
the number of particles and use these to compare the exact with 
the approximate dynamics. 



1. Introduction 

The purpose of our present work is to investigate certain aspects of 
the evolution of a large number of indistinguishable Quantum particles 
(Bosons) under binary interactions. If we call if)(t, x±, x-i . . . xn) the 
wavefunction describing the N particles with Xj G M 3 the coordinates 
for j = 1, 2 ... N, then tp satisfies an evolution equation of the form 

(l/W = {Hi - ArV) (1) 

where H 1 is a sum of the from Y^j=i ■ The term V models two body 
interactions of the following general type 

V := (1/2) N 3p v(N p { Xj - x k )) ; < (3 < 1 (2) 

where v is positive, even, bounded and integrable. 

In the equation (CTJ above we consider non-relativistic particles and 
set h = 2m = 1 for simplicity. Here and for the rest of this pa- 
per we denote vn '■= N^^v^N 13 ■ ). The fact that we consider Bosons 
means that the wavefunction is invariant under all permutations of the 
indices j = 1,2 ... N and one would like to solve the evolution equa- 
tion under some initial condition at, say t = 0, ^(0, Xx, x 2 . . . x n ) : = 

l 
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ijjo(xi,X2 ■ ■ . xjv). Presently we are interested in the evolution of fac- 
torized (or approximately factorized) initial data i.e. we would like to 
consider special initial data of the form 



The evolution of ([Q) with initial data ([3]) is quite complicated for iV 
large and one would like to have an effective approximate description 
of the evolution. The motivation for this type of problem comes from 
Bose-Einstein condensation where one considers a large number of iden- 
tical (indistinguishable) particles in a trap. Einstein following ideas of 
Bose, observed that nonineracting particles in a box undergo a phase 
transition at a critical temperature proportional to density 2//3 , so that 
below this temperature a macroscopic number of particles occupy the 
ground state, furthermore at zero temperature all particles condense 
to the ground state. It is natural and more realistic to consider inter- 
acting particles. Following ideas of Landau a heuristic theory based on 
the idea of the mean-field approximation was developed by Gross and 
Pitaevski [32], [22]. On a more fundamental level, the problem of a 
weakly interacting Quantum gas was taken up in the pioneering work 
of Lee, Huang and Yang as well as Dyson) [27] , [10] . 

More recent theoretical developments are due to Lieb, Solovej, Yng- 
vanson, Seiringer et. al. see [2H] and references therein. In particular, 
Theorems 6.1, 7.1 in [29], as well as Appendix C in [TS] strongly suggest 
that that the ground state is well approximated by a tensor product 
fl3]) where 0o describes a mean field approximation. 

Let us point out that we can in fact treat more general initial con- 
ditions corresponding to the iVth component of e ^N A (<t>) e B ( k ) see 
section (T5]). 

Experimental confirmation of Bose-Einstein condensates was finally 
achieved [T], using alkali atoms. The reason for the use of alkali atoms 
is the fact that they contain a single valence electron in the outermost 
s-orbital (for example 5-s for Rubidium). The other contributions to 
the total spin comes from the nuclear spin. If the nuclear isotope is one 
with odd number of protons and neutrons it will have a net half integer 
spin. For example Rubidium 87 has S = 3/2 from the nucleus. The 
total spin takes the values S = 1 or S = 2. If we prepare the sample so 
that only one of these states is present then this will be a gas of identical 
Bosons. If two different states are present then we should consider it 
as a mixture of two different gasses. Since alkali atoms are complicated 
composite particles their interactions are not known explicitly, which 
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means that the potential v in (j2J) is not explicitly known, moreover one 
can treat the atoms as Bosons only for a sufficiently dilute gas. At 
shorter distances the internal structure of the atoms should be taken 
into account. Here we consider a sufficiently dilute Boson gas and we 
make the reasonable assumption that the interactions are repulsive i.e. 
v > and that they are short range in the sense that J v(x)dx < oo. 
It is clear from the above comments that one should consider particles 
with spin. The present framework can be generalized in this case in 
a straightforward manner, however in the name of simplicity we forgo 
this generalization. 

Let us comment on the scaling present in the form of binary interac- 
tions. The parameter describes the strength of particle interactions. 
It is a reasonable (but not obvious) idea to assume that the evolution 
of is approximated by a tensor product i.e. 

N 

and the issues are, first to explain the nature of the approximation 
described in P| and second to derive an evolution equation for 4>(t,x) 
consistent with the dynamics of (pQ). On the second question, the gen- 
eral idea is that the evolution of the mean field <f>(t, x) satisfies an 
equation of the form 

id t <f> = A0-^(|0| 2 )0, (5) 
and the nonlinear term gp(\cf)\ 2 ) depends on (3 in the following manner, 

So(|0D = f dy{v(x-y)\<t>(t,y)\ 2 } 

9M*) = (Jv(y)dy\\ ( j ) (t : x)\ 2 ; < < 1 

g 1 (\ ( j ) \ 2 ) = 87ra\ ( j ) (t,x)\\ 

where a appearing in g± is the scattering length corresponding to the 
potential v. In the case (3 = one obtains a Hartree type evolution for 
the mean field and considerations similar to our present work where 
taken up in [20,], [21]. The case j3 = 1 is probably the most interesting. 
In this case the scaling is critical in the sense that particles develop 
short scale correlations which in the limit iV — > oo lead to the appear- 
ance of the scattering length in the equation. A heuristic argument for 
this is well known in the Physics community, however the explanation 
on how the scattering length emerges from the N body dynamics was 
recently given in the work of Erdos, Schlein and Yau [16].|15j. 
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Our aim is to introduce pair excitations as a correction to the mean 
field approximation. This goal is achieved by introducing a kernel 
k(t,x,y) which describes pair excitations and one would like to derive 
an evolution equation for k consistent with the N body dynamics, 
which means that we should be able to obtain estimates comparing 
the exact with the approximate dynamics. The general idea of the 
approximation can be described in the following manner. Two particles 
leave the condensate and form a pair v^{x\ — £2)0(^1)0(^2) which in 
turn drives the evolution of pair interactions. It turns out that a natural 
way to introduce pair excitations as a correction to the mean field is via 
a Fock space formalism which we will outline in the next section. Let 
us comment here on the nature of our approximation. The mean field 
approximation (j4]) is a simple description of the iV-body wavefunction, 
however the nature of the approximation is quite involved and uses the 
BBGKY hierarchy and its limit as N — > oo as shown by Elgart, Erdos, 
Schlein and Yau [TT1 112 H3 13 E51 US] • See the approach of [25], [26], 
[5] based on space-time estimates. We also mention the related case of 
3 body interactions [6] , and switchable quadratic traps [S] . 

Moreover the approximation does not track the exact dynamics, 
rather its true usefullness lies in the fact that it can (approximately) 
track observables. In contrast our approximation is more complicated 
but it tracks the exact dynamics in Fock space norm. As a matter of 
fact a heuristic explanation of our approximation runs as follows: The 
iV-body wave function consists of three parts, particles that live in the 
condensate, bound pairs and particles that decayed after forming pairs. 
Controlling the number of particles that formed pairs leads to another 
justification of the mean field approximation. We will not pursue this 
line of inquiry here, however the approximation can be readily used to 
estimate observables. There are two main points in our present work. 
First we have a new transparent derivation of the evolution equation 
of pair excitations, indeed we derive a new system of linear equations. 
Second we obtain apriory estimates for the pair excitations kernel which 
are independent of iV and this, in turn, allows us to estimate the dif- 
ference between the exact and approximate solutions provided that /3 
is sufficiently small (/3 < |). 

Our work was inspired by [35] as well as [37]. Previous works directly 
related to the present are [H] and [24] . See also [2] and [3Tj . 

See Theorem (12. 3p below for the precise statement of our main result. 

The paper is organized as follows. In section 1 we develop the Fock 
space formalism which is necessary for our computations and derive 
the evolution equations for the pair excitation kernel. In section 2 
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and 3 we we derive the apriory estimates for the mean-field and for 
the pair excitation kernel. In section 4 we show how this information 
can be implemented in order to compare the exact solution to our 
approximation. 

2. FOCK SPACE FORMALISM AND THE NEW DERIVATION 

In this section we introduce the Fock space formalism and the Hamil- 
tonian evolution in symmetric Fock space. F is a Hilbert space consist- 
ing of vectors of the form 

= (ipo , Vi(zi) , ^2(^1,^2) , • • • ) 

where ipo e C and ipk are symmetric L 2 functions. The norm of such a 
vector is, 

00 

11 wir=w>=w 2 +E 

71=1 

Thus F is a direct sum of sectors F„ of the form, 

00 

F = ^F n ; ¥ n :=Ll(R 3n ) 

n=0 

with Fo = C and L 2 S (M. 3 ) denoting the subspace of symmetric functions. 
In the Fock space F we introduce creation and anihilation distribution 
valued operators denoted by a* and a x respectively which act on sectors 
F n _x and F„ +1 in the following manner, 

1 " 

a*(V>n-l) :=^^^(x- X^n-iixi, Xj-i, Xj+i, ...,x n ) 

, X\ , . . . , x n ) 

with [x] indicating that the variable x is frozen. In addition a x kills F 
i.e. a x (ipo) = 0. The vacuum state will play an important role later 
and we define it as follows 

|0> := (1,0,0...) 

so that a x \0) = 0. One can easily check that [a^a*] = 5(x — y) and 
since the creation and anihilation operators are distribution valued we 
can form operators that act on F by introducing a field, say (f)(x), and 
form 

:= J dx {4>(x)a x } and a*^ :— J dx {(j){x)a* x } 

where by convention we associate a with and a* with 0. These 
operators are well defined, unbounded, on F provided that is square 
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integrable. The creation and anihilation operators provide a way to 
introduce coherent states in F in the following manner, first define 



A(4>) := J dx {(j)(x)a x — 0(x)a*} 

and then introduce iV-particle coherent states as 

|^(0)) := e -^W|0) . (6) 



It is easy to check that 

e -^W|0) = ^. . . Cnfl^xj) ...) with c n = (e~ N N n /n\) 



1/2 



In particular, by Stirling's formula, the main term that we are inter- 
ested in has the coefficient 

c N « (2ttN)- 1/4 (7) 

Thus a coherent state introduces a tensor product in the sector Fjv, 
hence we can use such states as a mean field approximation to the 
Hamiltonian evolution in Fock space, see fl4]). 

The Fock Hamiltonian (acting on Fock space vectors) is 

U:=U X - N~ X V where, (8a) 

"Hi := J dxdy {A x 5(x — y)a* c a y } and (8b) 

V := ~ J dxdy {v N (x - y)a* x a*a x a x } , (8c) 

where we set 

v N (x-y) :=N^v(N^\x-y\) , (9) 
and the evolution in Fock space is described by the equation, 

-d t \^) = U\tlj) (10) 

% 

which has the formal solution 

\^{t)) = e m \%) . (11) 

Notice that % preserves the sectors F n and that H agrees with the 
classical Hamiltonian (pQ) on Fjy. However in this framework we allow 
any number of particles to evolve and one is interested, in particular, 
in the evolution on the sector Fjv. 
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Our goal is to approximate \ip(t)) in (ITTj) and for this purpose we 
introduce two fields 4>(t,x) and k(t,x,y) and the associated operators, 

A{4>) :— J dx {(p(t, x)a x — (p(t, x)o* x } (12a) 

B{k) := dxdy {k(t, x,y)a x a y — k(t, x,y)a x a*} . (12b) 



The coherent initial data are introduced via J^o) — e - ^' 4 ^ -' |0) 
which means that the initial data are a tensor product on Fjv as desired, 
see ([3]). Our approximation scheme is 

|Vw> := e-^We-sVe™*®^) (13) 

with x(t) a phase factor, and we plan to show that |V ; (^)) ~ |VVpr(^))- 
The issue for us is to determine the dynamics of the fields <fi and k. 
It turns out the the evolution of k is described via a set of new fields, 

sh(/c) := k + ^k o ko k + . . . , (14a) 

ch(A;) := S(x - y) + ^k o k + . . . , (14b) 

where o indicates composition, namely k o / stands for the product, 

(k o l)( Xl x 2 ) := / dy {k(x ll y)l(y,x 2 )} . 



A crucial property of the above multiplication is that it is not com- 
mutative i.e. k o I I o k. In order to describe the evolution we need 



g N (t,x,y) := -A x 8(x-y) + (v N * \<f)\ 2 )(t, x)5{x - y) 

+ v N (x - y)4>(t,x)(f>(t,y) (15a) 
m N (x,y) := -v N (x - y)(f)(x)(p(y) . (15b) 

Using y jy we can construct two operators as follows: For a function 
s(t,x,y) symmetric in (x,y) and a function p(t,x,y) conjugate sym- 
metric in (x, y) i.e. p T = p, we define 

S(s) := -s t + gjr o s + s o g N (16a) 

i 

W(p)--=~Pt + [gJf,p] ■ (16b) 



s 
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The dynamics of the fields are determined via, 

-d t (j)- A0+ (v N * |0| 2 )0 = O (17a) 



S (sh(2Jfe)) = m N o ch(2k) + ch(2k) o m N (17b) 
W (ch(2fc)) = m N o sh(2fc) - sh(2A;) o . (17c) 

Remark 2.1. It is clear that ch(2/c) and sh(2/c) are not independent of 
each other, thus we can ignore the third equation, however in the form 
stated above the equations are readily amenable to the derivation of 
apriory estimates. The equation for <fr is of Hartree type and its formal 
limit as iV — > oo is NLS. 

The theorem concerning the evolution of the mean field <fi and the 
pair excitation kernel k reads as follows. 

Theorem 2.2. Suppose that < j3 < 1 in (Q. Given initial data 
</>(0, x) := (po(x) G W k ' 1 (k derivatives in L , with k sufficiently large) 
and k(0,x,y) : = the system ( I17ap(117b|) ( ll7cp has global solutions in 



time which satisfy the apriori estimates, 

I //-i : n < C s (18a) 



C 

»(R8) + ||<9t0(£)IU°°(R3) < Tjj/a ( 18b ) 

sh(2fc)(t)|| x2(R 6 ) + ||ch(2A;)(t)-5|| i2( M6) <Clog(l + t) . (18c) 



The main difficulty in obtaining the estimates in the theorem above is 
the fact that defined in (jUJ) has a formal limit vn(x — y) — > c5(x — y) 
which means that has a limit which is not square integrable, as 
a matter of fact it does not belong to any L p for p > 1. In view of 
the theorem above, we can compare the exact with the approximate 
evolutions and the result is the following theorem. 

Theorem 2.3. Suppose that \ijj(t)) is the solution of (flOl) with initial 

data |V>o) : = e - ^" 4 ^ ) |0) and \il> a ppr{t)) is the approximation in ffTB"]) 
where the evolution of the fields <fi and k is determined from theorem 
(I2.2p . Under these conditions the following estimate holds, 

provided < (3 < | . This is a meaningful approximation of the Nth 
component of \^ provided^ < j3 < |, because of formula (I7j). A slightly 
more precise form of the estimate could be obtained by integrating the 
right hand side of the inequalities in Proposition (16. ip . 
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Remark 2.4. The real phase factor x is described via x(t) '■= /* ^i{A i o(^i) + 
iV -1 ^!^)} where 

M*) = -J / ^{^(a;-y)|0(t,x)| 2 |^Z/)| 2 } (20) 



2 _ 

and Hi is a complicated integral given in 

Proof. Here is an outline of the proof of this theorem. In order to 
relate the exact with the approximate solution we introduce the reduced 
dynamics 

|Vw(t)) := e VNA(t) e B(t )e itH e -VNA(0)^ (21) 

i.e. we follow the exact dynamics for time t and then go back following 
the approximate evolution. Notice that |VVed(0)) = |o) and if our 
approximation was following the exact evolution we would have that 
\ipred(t)) = |0). Thus our goal is to estimate the deviation of the 
evolution from the vacuum state. It is straightforward to compute the 
evolution of |VVe<f) an d it is 

-d t \lpred) = Ured\Aed) (22) 

where the (self-adjoint) reduced Hamiltonian is, 
1 



n red :=-(d t e B )e- e 

i 



The main idea is that the evolution of the fields <p and k is chosen so 
that the reduced Hamiltonian looks like 

U red = N/i(t) + J dxdy{L(t,x,y)a* x a y } - N- 1/2 S(t) 

where S(t) is an error term containing polynomials in (a, a*) up to 
degree four, and L is some self-adjoint expression which is irrelevant 
for the rest of the argument. 
Next consider 



:= e~ iNx(t) {\ij red )) - |0> where X (t) := / v(h)dti 
where we called \i := /x + N' 1 ^!- Thus 

-d t - -H re , + iV/i(t)^ e- w *W|Vv«i) = 0. 
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and therefore 

-A - Hred + Nfl{t)j |^> = N-V 2 e(t)\0) 

The equation above has a forcing term namely N~ 1 / 2 E(t) |0) and a 

standard energy estimate together with the fact that and e B are 

unitary, gives 

|| - |Vw(*))|| F 

= \\\m)\\ ¥ < N^ 2 [ t dtl {\\S(t l )\0)\\ w } . (24) 

JO 

The proof will be complete if we estimate the right hand side in the 
above inequality. Notice that £|o) has entry only in Fock sectors Fj 
for j = 1, 2, 3, 4 and in order to estimate it we need the lemma below 

□ 

Lemma 2.5. The error term is described as follows, 

E:=e B (Va + N-WVt) e~ B 

where V3 and V4 are cubic and quartic polynomials in (a, a*) respec- 
tively. Moreover the following estimate holds ifO < j5 < | ; 

||£(t)|0)|| F <CiV 3/3 / 2 log 4 (l + t) . (25) 

A more precise estimate is given in Proposition (16. ip . 

Remark 2.6. The polynomials P3 and P4 appearing in the error term 
are given by the expressions, 



-p 3 : = J dxdy {v N (x - y) (<j)(y)a* x a* y a x + 4>(y)a* x a x a y ) } 
P 4 : = (l/2) / dxdy {v N (x — y)a* x a*a x a y } 



as we will see shortly. 

The rest of this section is devoted to the derivation of (I17al) (117bl) ()17cj> . 
We have to compute H re d above, see (1231) . and for this task there are 
two crucial ingredients. They are based on the formal identities below 
for any two operators, say A and H, 

00 1 

e A He' A = -7 (ad^)" (H) (26a) 

n=0 n ' 

00 

(dte A )e- A = J2M^r 1 (A t ) (26b) 



ni 

n=l 
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where ad^("H) := LA, "Hi. They indicate that we have to compute 
repeated commutators of various operators. The series defining the 
exponentials in (I26aj) . ( I26bl) converge absolutely on the dense subset of 
vectors with finitely many nonzero entries provided that A = A(<f)) is a 
polynomial of degree one with <p 6 L 2 or A = B(k) is second order with 
1 1 k | \l2 small. If B is skew- Hermit ian, e B extends as a unitary operator 
for all k e L 2 . This construction is closely related to the Segal-Shale- 
Weil representation, as explained in [21], [17], [36], and our appendix 
([7j). This calculation was also used in our previous papers [201 f2T| . 

The first observation is the fact that since A(4>) is a degree one poly- 
nomial, if we denote by V n a homogeneous polynomial of degree n then 
commuting with A produces [A, V n ~\ = V n -i i.e. a homogeneous poly- 
nomial of degree n — 1. This in turn implies that repeated commutators 
produce a finite series in (I26ap . (126bj) which can be computed explicitly 
after some tedious but straightforward calculations. 

The second observation is that in (I26al) . (126bl) when we replace A 
with B we obtain infinite series with a certain periodicity which al- 
lows for explicit summation. This can be expressed via a Lie algebra 
isomorphism. For symplectic matrices of the blocked form 

'd(x,y) l(x,y) 
k(x,y) -d(y,x) 

where d, k and I are kernels in L 2 , and k and I are symmetric in (x,y), 
we define the map from L to quadratic polynomials is (a, a*) in the 
following manner, 



L :-- 



(27) 



The crucial property of this map is the Lie algebra isomorphism 

[^(LO,!^)] =X([L 1 ,L 2 ]) (28) 

thus any computation that involves commutations can be performed in 
the realm of symplectic matrices and then transfered to polynomials in 
(a, a*). In particular if we call I(H ) = 7-L for a quadratic Hamiltonian 
and I(K) = B then we have the two formulas below, 

e B Ue- B = l(e K He- K ) (29a) 

{d t e B )e- B = l((d t e K )e- K ) . (29b) 

Actually, to avoid the infinite trace in (I29ap . we write 

e B Ue~ B = U+[e B ,U}e- B 

= n + l([e K ,H]e- K ) 
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As a matter of fact if we define the following quardatic expressions, 

^xy - = a x a y j T^xy - = a x a y 
Q-xy := a x a y ) Qxy := a x a y 

then we can write, 

X(L) = ~ J dxdy{d(x,y)Vl y + d(y,x)V X y + k(x,y)Q* xy -l(x,y)Q X y} 
Remark 2.7. Notice that D\ = D yx + 8{x — y) thus we can write 
X(L)=-i J dxdy {d(x,y)V yx + d(y,x)V xy + k(x,y)Q* xy - l(x,y)Q xy ] 

— - J dx{d(x, x)} . 
In our present formalism if we define the matrix 



k o 



then we have that T(K) = B, see the expression in (112bl) . The expo- 
nential of K can be computed, 



6 = vsh(fc) ch^yj where ' 

sh(fc) := k + —k o k o k + . . . , 

ch(fc) := 5(x — ?/) + ^y/c o k + . . . , 

and o indicates composition. For completeness and for the convenience 
of the reader we include in the appendix the derivation of ( BHD , see a ^ so 

Let us now proceed with the calculations. First look at the expression 
inside the parentheses in the reduced Hamiltonian f[2"3"j) . It is straight- 
forward after repeated (but finite) commutations with A to come up 
with the expression below (see section 3 of [20]). 

-(d t e VNA )e- VNA + e^ A He-^ A 

i 

= Nfi + N X ^V X +V 2 + iV" 1/2 P 3 + N~ l V± (30) 

where V n indicate polynomials of degree n to be given explicitly below. 
The first term /xq is a scalar which can be absorbed in the evolution as 
an extra phase factor. It is given by the commutators, 

i [A, d t A] + \ [A, [A, H}] - i [A, [A, [A, [A, V 
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which reduce to the expression below, 

Hq := J dx |^t(#« - <Mh) - |V0| 
1 

~ 2 

The first degree polynomial V\ arise from the commutators, 



dxdy{v N {x-y)\<\>{x)\ 2 \<\>{y)\ 2 } . (31) 



+ [A, H]-j { [A, [A, [A, V 
and it can be expressed as follows, 

V 1 = J dx {h(t,x)a* x + h(t,x)a x } (32) 

where h := —(l/i)d t <p + A<p— (v^* |0| 2 )0- The second degree polyno- 
mial consists of the terms 

U x -\ [A, [A, V]] 

and can expressed 

V2 = \j dxd y^~9N{t,x,y)V yx - g N (t,y,x)V Xjy } 

+ \ J dxd y { m N(t,x 1 y)Q xy + m N (t,x 1 y)Q x y } (33) 

where g n and rriN are given by, see ( 115ajl . ( j!5bl) 

g N (t,x,y) := -AJ(x - y) + (fjy * \<p\ 2 ){t,x)5(x - y) 

+ v N (x - y)0(t,x)(j)(t,y) 

m N (x, y) := -v N (x - y)<f>(x)<p(y) . 

It is clear that (?at and above depend on the number of particles N. 
Subsequently, for simplicity, we will suppres this subscript and recall it 
only when it is relevant in an argument. Let us define the two operators 
below 

^ G := \ J dxdy {~ gN & x ' y^v* ~ gN & y-i x^X'-y} ( 34a ) 

M -- = \J dxd y { Tn N(t,x,y)Q xy + m N (t,x,y)Ql y } 

so that we can write P 2 = Hg + M.. The relevance of this splitting 
will become clear shortly. The third and fourth degree polynomiasl 
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arise from the commutators \A, V] and V respectively and are given 
below 



V 3 ■= J dxdy {v N (x - y) (<f)(y)a* x a* y a x + <f)(y)a* x a x a y ) } (35a) 
V4 := (1/2) / dxdy {vj^(x — y)a*a*a x a y } . (35b) 



The mean field approximation emerges from the first degree polynomial 
V\. Since /io can be absorbed into the evolution it is reasonable to pick 
the field so that h{4>) = 0. This leads to the evolution 

-d t 4>- A<f)+ (v N * |0| 2 U = O (36) 
1 

which is of Hartree type. The formal limit of the equation above is 
the cubic NLS where the constant in front of the nonlinear term is the 
integral of the potential v. If (f> satisfies ( 136]) then /x reduces to 

fio = ~ J dxdy{v N (x-y)\(f>(t,x)\ 2 \(f)(t,y)\ 2 } . (37) 

Now we can compute the reduced Hamiltonian in fl23l) using the 
splitting in (134aj) . f!34bj) . First let us first give a name to 

£ : = e B (V 3 + N-^ 2 V 4 ) e~ B (38) 

which will be treated later as an error term. Now we can write, see 



n red = - (d t e B ) e~ B + n G + [e B , U G \e~ B + e B I(M)e~ B 
- e B (N-^Vs + N- l V A ) e- B 

= U G +X ((1/i) {d t e K ) e- K + [e K , G]e~ K + e K M e - K ) - N~^ 2 S 
= n G +I(R) - N- 1/2 S , (39) 

where R is defined to be the expression, 

R := (1/i) (d t e K ) e- K + [e K , G]e~ K + e K Me~ K . 
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For the convenience of the reader, let us recall our set up, 



K-(° k ) and e K -( cHk) sh( ^ 
K -{k Oj and 6 -{sh(k) ch(k)) 

sh(k) := k + —k o k o k + . . . , 

ch(A;) := 5(x — y) + -^k o k + . . . , 

g(t, x, y) := -AJ(x - y) + (v N * |0| 2 )(t, x)5(x - y) 

+ v N (x - y)4>(t,x)(j)(t,y) , 
m(x,y) := — Vjy(x — y)(j)(x)(f)(y) where vn(x) = N 3 ^v(N^x) 

n °t) and M := ( ° f 
— g / V— m 



N u : = ^ (this corresponds to the Number operator) 

S(s) := -s t + g T o s + s o g and W(p) := + [g T ,p] . 

Thus S describes a Shrodinger type evolution, while W is a Wigner 
type operator. These operators will emerge shortly. Recall the formula 
f )39|) that we derived earlier for the reduced Hamiltonian 

'Hre d = U G + Z{R)-N- 1 l 2 E , 

where %q has only a* a terms (which annihilate the vacuum) and R 
can be computed explicitly. In fact we have, 

R = 



1 fdx(k) t sh(fc) A 
z Vsh(A;) t ch(&)J 

/ [ch(fc), <?] -shorn — sh(/c) o g T — g o sh(/c) + chom 

\sh(fc) o g + g T o sh(k) -chom — [ch(fc), g T ] + sh o fn 



o , ch(fc) -sh(fc) 



-sh(fc) ch(fc) 

(matrix product, where kernel products mean compositions) 

The condition that we would like to impose is that R is block diagonal 
so that T(R) contains only terms of the form aa* and a* a so that, apart 
from a trace when we commute a with a*, we obtain an operator which 
annihilates the vacuum state. The remaining trace can be absorbed in 
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the evolution as a phase factor. Thus our requirement is 



d 

d~t 



e K ) e~ K + [e K , G]e~ K + e K Me~ K is block diagonal. (40) 



-if 



-K 



K 



We proceed to show this equivalent to equations (117bj) . (I17cl) . Let us 
make the elementary observations 







dt 



k d _ 
e o — o e 
dt 



K 



G 



dt 

[e K ,G]e- K = e K Ge~ K 

so removing the part of ( 1401) that is diagonal already we have the 
equivalent formulation of ( 1401 

1 d 



K 



i dt 



+ G + M ) e K is block diagonal. 



(41) 



Now we make the observation that a matrix is block-diagonal if and 
only if it commutes with the number operator matrix N u , as well as 
(for arbitrary matrices A and B) we have [e K Ae~ K , B] = if and only 
if [A, e~ K Be K ] = 0, so our equation (j4Tl) reads, 

1 d 



i dt 

A direct calculation gives 



G + M 



N u e> 







(42) 



-K 



K 



ch(2fc) sh(2k) \ 
-sh(2Jfe) -ch(2fc) J 
after which is is straightforward to compute 



ld_ 

"Tat 

and simlarly, 



+ G , e~ K N u e 



W ch(2fc) S((sh(2fc)) 



S (sh(2fc)) W ch(2A;) 



[M, e~ K N u e K ] 



-m o sh(2A;) + sh(2A;) o m —m o ch(2A;) — ch(2fc) o m 
—m o ch(2A;) — ch(2fc) o m —m o sh(2fc) + sh(2A;) o fn 



Finally combining the two formulas above we obtain, see (j4"2"|) . the linear 
pair of equations below 

(43a) 
(43b) 



S (sh(2fc)) = m o ch(2fc) + ch(2A;) o m 

W (di(2k)} =mo sh(2k) - sh(2k) o fn . 

This completes the derivation of the evolution equations for the pair 
excitations and the mean field, namely (T43aj) . (l43b|) . together with ([36D 
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describe the evolution of and k and are the equations in ( ]17a|) . (117bj) 
and (I17cp . In particular, we have proved that in that if 0, k satisfy 
these equations, then the energy estimate ff24"j) holds. 

3. Estimates for the solution to the Hartree equation 

This section adapts classical results for NLS due to Lin and Strauss 
[3D] . Ginibre and Velo [T9], Bourgain [3], as well as Colliander, Keel, 
Staffilani, Takaoka and Tao jl] to the Hartree equation. Assume 

I^0_A0+(^*|0| 2 )0 = O (44) 
i at 

0(0, •) = 00 • 

We recall the relevant conserved quantities, following the notation 

P'= (V2)|0| 2 ; 

Pj := (1/2^ (0V.-0 - 0V,0) ; po = (l/2i) (0^0 - 0^0) ; 
o-jk ■= Vj0Vfc0 + V fc 0Vj0 ; a 0j = V,-0<9t0 + d t 4>Vj(j) 

A:=-S(09 t 0) + |V0| 2 + i(t;*|0r)l0| 2 

= ^(A|0I 2 -(^*I0I 2 )I0I 2 ); 

e: = \Vd>\ 2 + ^(v*\d>\ 2 )\(f>\ 2 . 
The associated conservation laws are 

dtp ~ Vy' = , (45a) 

d t Pj - V, | a/ - ~5/aJ + = , (45b) 

<9 4 e - V i( r Q j + l = . (45c) 

These laws express the conservation of mass, momentum and energy, 
respectively, where the vector (lj,lo) is 

lj := 2 ((tijv * Pj)p - (% * p)Pj) , /o := 2 ((t>jv * Po)p - (vn * p)po) 

We adapt the well-known method of interaction Morawetz estimates, 
due to Colliander, Keel, Staffilani, Takaoka and Tao, outlined in [I]. 
Start with 

Q{t) - / fojP i ( t i x )p{t, y) + P(*> ^) VjP J (t, x)) |x - y\dxdy . 
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Using ( 145]) we get 

Q(t) = 2 J V/(t,x)V fc /(t,t/)|i - y\dxdy 

+ J (vJvJa/^x)-^^^^)^-^^)^ p(t,y) 

+ p(t, x)Vj ( [V k jcx/(i, y) - ^/A(t, y) j - y)^ ) J I* ~ 3/1^2/ 

> ^ f _ ^ A (*> a; )/ ? ( t ' ~ ^(*> S/)J ^l 35 ~ 2/1^^2/ (main term) 

((Vjlj(t,x))p(t,y) + p(t,x)(Vjlj(t,y))) \x — y\dxdy (error term.) 

We have used the fact which we recall for the reader's convenience (see 
[4J), that 

(Vj Vfca) (x - y) (-2pj(t, x)p k (t, y) + a/(t, x)p(t, y) + a/(t, y)p(t, x)) 



(Vj-Vfco) (z - y) + 4>j{x)(t){y)) (<t>{x)4>j{y) + </>j(x)<j>(y)) 



+ {<f>{x)<t>j{y) - <j>j{x)<t>(y)) (<f>(x)<l> } (y) - ^(x)0(y))J > 

where a(x) = As it is well known, 

(main term) > c||0(t, -)IIl 4 w ith c > 

while the error term is in fact 0, by inspection. Since Q(t) is bounded 
uniformly in time by ||0o||#-i> we have shown the following proposition. 

Proposition 3.1. Let <fi be a solution to the Hartree equation 
There exists C depending only on H^oIIh 1 such that 



|L 4 ([0,oo)xIR 3 ) < C 

and, as an immediate consequence of conservation of energy, 

||0||l8[o,oo)l4(k3) < C . (46) 

Remark 3.2. It was shown by Bourgain in [3] that if is a solution to 
cubic NLS, then there exists C s depending only on ||0o||i? s such that 

\m,-)\\ H s<c s vt . 

Using the above Morawetz estimate (which was not yet discovered 
when Bourgain did this work), we can easily prove the same for our 
Hartree equation. 
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Proposition 3.3. Let <fi be a solution to the equation ( 1441) . There 
exists C s depending only on \\4>q\\h s such that such that 

\m,-)\\B*<c s 

uniformly in time. 

Proof. Split [0, oo) into finitely many intervals Ik where 

||0||L8(/ fc )L4(K3) < e 

where e is to be prescribed later. Differentiating (|44|) 

j J^V - AD S <P = -D s ((v N * |0| 2 )0) (47) 
where 

D s ((vn * |0| 2 )0) = (vn * |0| 2 ) D s (f) + similar and easier terms. 
For the first interval, ii, we get, using the L 8//3 L 4 Strichartz estimate, 
\\ D '<P\\ifi/*{i 1 )L*(a*) < C||0o||^ + C\\ (v N * |0| 2 ) D s 4>\\ Lm{ i l)L i, 3m 

< Cl||0o||H= + C2 1 1 1 1 i 8 (/i)Z, 4 (M 3 ) 1 1 1 1 Z, 8 / 3 (Zi)Z,4(K.3) • 

At this stage, we pick e so that C 2 e 2 < | to conclude 

In turn, this allows us to control the inhomogeneity of (|4"Tj) 

|| (V N * |0| 2 ) ^^11^/5(^)^/3^3) 
^ I|0|Il8(7 1 )L4(] R 3)||-D S 0|| jL 8/3( /i ) L 4( R 3) 

<C||0o||^ 

and therefore 

-)\\h' < \\Mh° + C\\^l 8{Ii)L4m \\D s ^ L 8/3 (h)L 4 m 

< c\\Mhs 

for all t E I\. Repeating the process finitely many times, we are done. 

□ 

If we assume the data <fto and sufficiently many derivatives are not 
only in L? but also in L 1 , we can also get decay. 
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Corollary 3.4. Let <p be a solution to (jUJ). There exists C depending 
only on ||0o||wm f or k sufficiently large such that 

||0(V)IU»<- (48a) 
and also 

IW(V)IU-<^ (48b) 

Proof. We solve f|44l) by Duhamel's formula and use the standard L°° L 1 
decay estimate for the linear equation. We notice the following trivial 
estimates, based on the fact that \\4>(t, -)\\h 2 is uniformly bounded: 

\\e^ A (v * |0| 2 0( S )) lUoo < C\\ (v * |0| 2 0( S )) Hj^ < C||0( S , 
We have 

||0(V)|U~< A||0o|Ui+ / u „L , J |0(3,OIU~tfa 



t 3/2 II WJII^ | t _ s | 3 /2 + 1 

or, denoting M(t) = (1 + t 3/2 )||0(t, -)IU<», We have 

This implies M(t) < 2C for t sufficiently large. A similar proof works 
for the time derivative. □ 

By interpolating with the L 2 uniform bound we get the next Corol- 
lary. 



Corollary 3.5. Let <p be a solution to fl44j) . There exists C depending 
only on H^ollvy*- 1 f or k sufficiently large such that 

||0(V)|| L 3 + ||d t 0(V)||L3<-^ ■ 

1 + 



4. Estimates for the pair excitations 

Define ch(2/c) := 5 + p 2 , sh(2A;) := s 2 , and also ch(fc) := 8 + pi, 
sh(2A;) := Si so that, see (143 aft . (143bj) become 

S (s 2 ) = 2m + m o p 2 + p 2 o m (49a) 

W (p 2 ) = m o s 2 — s 2 o fn (49b) 

S2 (0,-)=p 2 (0,-) = 
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The goal of this section is to prove the following theorem. 

Theorem 4.1. Assume <fto £ W k,x for k sufficiently large. The follow- 
ing estimates hold: 

\\s 2 (t, Oilmen..) + \\P2(t, -)\\vqp) < Clog(l + t) (50) 



where C depends on ||0(O, for some finite k. A similar result 

holds for the higher time derivatives, but we will not use it or prove it. 

An immediate corollary is of the above theorem is, 

Corollary 4.2. The following estimates hold: 

< Cl0g(l + t) 

lbi(V)IUw < cio g (i + t) 

\p-i_(x,x)\dx < CTog 2 (l + t) . 



Proof, (of corollary (14.21) ) Since sh(2A;) = 2sh(A;) o ch(k), we get 
\\si(t, -)IU 2 (K6) < 2 ll s 2(^, •) IU 2 (iR6) ||ch(A;) _1 || opera f 0r . 

< 2 N^aC*, -)IU 2 (» 6 ) • 

We also have pi(x,x) > 0, p\ opi(x,x) > 0, so taking traces in the 
relation 

Pl O p x + 2p x =8! OS! 

gives the other estimates. □ 

Before starting the proof of Theorem ( 14.1 ft , we need some prelimi- 
nary lemmas. 

Lemma 4.3. Recall m(t,x,y) = —vn(x — y)<p(t,x)<p(t,y). Then there 
exists C such that 

l ^^d^d V < CU(t,. )\\is (51a) 

(Ifr + \w) 

and also 

[ '^iY;^' 2 ^ < \m, whwm, oiiia • (sib) 



Similar estimates hold for higher time derivatives. 
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Proof. Write 

vn(x - y)(f>(t, x)(f)(t, y) = J 5{x-y- z)v N (z)(j)(t, x)<f>(t, y)dz . 

The Fourier transform of 5(x — y — z)<j)(t,x)<j)(t,y) is easily computed 
to be 

where we denote 6~(x) = Mx — z). Thus 



\fn£(t,£,Ti)\ = | J v N {z)e lz ^(j)4> z {t,i + ri)dz\ 
< \\ v n\\li / \v N (z)\\<f)(j) z (t,^ + r])\ 2 dz 



Thus, after a change of variables, the left hand side of (15 lap is domi- 
nated by 



\ v n(z)\ ^J z ^\^\ 2 d^dr]dz 
(If I + M ) 

< C J \v N {z)\ — didz 

<c ! \v N {z)\w\\U* = c\\4>\\» 



We have used the fact that 



z) Hi 2 



< C||00 2 || L 3/ 2 (by Hardy-Littlewood-Sobolev) 

< cupula ■ 

The proof of (151bj) is similar. □ 



Lemma 4.4. Let s° a be the solution to 

f~-A R 6^s° a {t,x,y) = 2m{t,x,y) (52) 
s o a (0,x,y) = 0. 
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Then 

ll s a(^ -)IU 2 (K 6 ) 

<c(ll0(o,-)lli3 + ||0(t,-)lli3+ J*U(s r )\\ L 4d s <P(s,.)\\ L3 ds 

< aiog(i+t) 

where C depends only on H^oIIwm- 

Proof. Solving by Duhamel's formula we get 

ll S a(^ -)IU 2 (K6) 

< C \\ / e is (l«l 2 +N 2 )m( S ,e,ry)rf S || i2(R6) 
Jo 

. „/,, m(0^) , H m(t,£,ri) \ 

- VWTW hHm "W+W iR6) ) ( } 
+ || ^^,w)M^) rfs | lL2(R6) . (54) 

io 161+ M 2 

The terms (1531) are estimated directly by Lemma (14. 3p . For the last 
term (1541) move the norm inside the integral and use Lemma (14. 3 j) as 
well as Corollary (13. 5p . □ 

Lemma 4.5. Let s a 6e i/ie solution to 

8(s a ) = 2m(t,x,y) (55) 
s a (0,x,y) = . 

T/ien 

0IU a c«i a ) - c^si 1 + t ) 

where C depends on \\<f>o\\wk,i- 

Proof. Let V be the "potential" part of S, so that 

1 d 

S = -—-A R6 + V. 

% at 

The operator V is bounded from L 2 to L 2 , with norm bounded by 
U(t, OHicc < igs (by Corollary (E3J)). Explicitly 

V(u)(t, x, y) = ((v N * \0\ 2 )(t, x) + (v N * \(f)\ 2 )(t, y)) u(t, x, y) 
+ J vn(x — z)(f)(t,x)(f)(t, z)u(z,y)dz + J u(x, z)vn(z — y)4>(t, z)<f>(t,y)dz . 
Write 

s = s° + s 1 
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where s° is as in the previous lemma, 
1 d 

x, y) - A R es°(t, x, y) = 2m . 
Then si satisfies the equation 

SW) = -v Mfr, ■)) ■ (56) 

Both s° and are zero initially, and the estimate is already known for 
s°. . We apply energy estimates to (1561) . Recall 



W(4 o si) = S(4) o a j - 4 o S(4) = -V {s° a (t, •)) o si + si o V W{t, •)) 
Taking traces, we get 

d y 2 1 

7^II S Jl 2 (IR6) ^ 2 \\ V ( S a(^ •)) lli 2 (lR 6 ) ll 5 alU 2 (M6) 

c 

< 1 +t3 ll g a(^-)l|L 2 (lR«)||-SalU2(R6) 

Clog(l + 1) 1( 

Integrating, we get the estimate 

\\sl(t,-)\\ L , m <C 

which implies the claim. □ 

We are ready for the proof of Theorem (14. ip . 

Proof. Write s 2 = s a + s e where S (s a ) = 2m, as in the previous lemma. 
The kernels s e and p 2 satisfy the following less singular system: 

S (s e ) = m o p 2 + p2 o m (57a) 

W (p 2 ) = (m o — s a o fnj + mos7 — s e om . (57b) 

Using lemma ( 14. 5 p to estimate \\s a (t, -)\\l 2 and defining 

M(t, •) = m o~s^ — s a ofn 



we have 



S (s e ) = m o p 2 + p 2 o m (58a) 
W (p 2 ) = M + mos7 — s e om (58b) 
where 

Clog(l + t) 



||M(t,-)||^(R6) < 



t 3 



since 



C 



\m o s\\ L 2 < C\\(j)\\ 2 Loa \\s\\ L 2 < ,J \s\\l2 

1 + i d 
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Using the general formulas 



W(s o s) = S(s) os-so S(s) 
W(p op) = W(p) op + po W(p) 
and taking traces we get 

c 

< Y^Tp (lb 2 ( £ ' •)IU 2 (K6)||Se(*, -)IU 2 (K6)) + <7||Af (*, -)IU 2 (K 6 )lb2(^-)IU 2 (R 6 ) 

which leads to the desired estimate. Explicitly, define 
E 2 (t) = \\s e (t,-)\\ 2 L2m + \\p 2 (t,-)\\ 2 L2 

then 

^m)<c( Kl ^E{t) + \\M{t 1 -)\\ L , 

thus E(t) is uniformly bounded. Also, HP2IU 2 stays uniformly bounded 
and the logarithmic growth of || S2 Hz, 2 can be traced back to ||Sq||l 2 from 
Lemma (jOj) . □ 

5. List of all the error terms 

Our purpose in the present section is to compute explicitly all the 
error terms and show how they can be estimated. Fortunately there 
are only a few terms for which one should be carefull, the rest being 
easier to handle. 

Let us recall here briefly our basic strategy, which is to define 

U ^ ._ e B(k(t)) e VNAm)) e m e -^NA(4>(0))^ 
and Compute V-red SUch that (l/i)d t \ljjred) = 'Hred^red) ■ 

n red = -| ( e w») e-w)) 

1 ut ' 

+ e B(k{t)) ^}_^_ e VNAm))J e -VNA(<t>(t)) +e VNA(<i>(t))ff e -VNAm))j e ~B{k{t)) 

= Nno(t) (59a) 
+ N l ' 2 e B ^ J dx {h((f>(t, x))a* x + h((f)(t, x))} e - B< - m) (59b) 
+ \d ^ mt)) j e _ mt) + eB{m ( _ 1 ^ [A v]]] \ e _ mt)) (59c) 

_ N -i/2 e BW)) f[ A) V ] + N^vj e- B ^ . (59d) 



26 M. GRILLAKIS AND M. MACHEDON 

The function fio(t) and h((f)(t,x)) appearing in (I59al) .( l59b1) are given 
below, 

:= J dx - <Mh) - | V ^| 2 } - \ j dxdy{v N (x - y)\<p(x)\ 2 \(p(y)\ 2 } 

h((f)(t, x)) := --d t <P + A0 - (v N * |0| 2 )0 . 



As we know, (I59bj) is set to be zero due to the Hartree equation for 



and (I59cp which contains terms of order 0(1) becomes block diagonal 
by an appropriate choice for the evolution of k. We can compute the 
0(1) term in flTO . 

where 



/ := (S(sh(fc)) - ch(k) o m) o ch(fc) - (W(ch(fc)) + sh(fc)om) o sh(fc) 
w := (W(ch(/c)) + sh(k) o m) o ch(fc) — (S(sh(A;)) — ch(k) o m) o sh(k) . 

The evolution implies that / = thus 

'w T 



pED = u G + i 



-w 



Multiplying / on the right with sh(k) and using the identities ch(k) 
sh(k) = sh(k) o ch(A;) and sh(fc) o sh(k) = (ch(A;)) 2 — 1 we discover 



/ o sh(k) + w o ch(k) = W(ch(/c)) + sh(fc) o m . 

Using the formula (see [21] , lemma 3.1, stated in slightly different no- 
tation) 



-i 



W(ch(fc)) o ch(fc) +ch(A;) oW(ch(fc)) 



= ch(k) o mo sh(k) — sh(fc) o m o ch(k) 

we get 



w(y,x) = W(ch(A;)) o (ch(A;)) 1 + sh(k) o m o (ch(/c)) ' 

^(dT(fe))^omoI(fc) + s li(it)omo (dn^))' 1 ) + ^ W(chp-)), (chfT)) -1 
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If we write w(y,x) for the above we have the quadratic term and a 
trace 



dSHcj) =u G dxdy i w ( x i y) v yx + w (y> x ) v x y } (60) 

+ -Tr^(ch(A;)) 1 o m o sh(k) + sh(k) o m o (ch(/c)) 1 



The error terms, which are order 0(iV 1//2 ) or higher, are all in the 
third line (159dj) i.e. what we called £, 



S = e B ([A,V] + N- 1 / 2 V) e- B . 



Recall that 



[A,V] = J dx 1 dx 2 {v N {xi -x 2 ) (<fr(x 2 )a* Xi Q xlX2 + <fr(x 2 )Q* XlX2 a Xl ) } 



(61a) 



V = i y dxidx 2 - a^Q^Q^a} (61b) 



and the transformation of (a, a*) is 



e B a x e B = J (ch(k)(y, x)a y + sh(k)(y, x)a*) dy := b x (62a) 
J (sh(k)(y,x)a y + ch(k)(y,x)a^j dy := b* x . (62b) 



e B a*e- B 



One can see that [b x ,b*] = S(x — y). For the rest of this section, 
the argument of the "trigonometric" functions is always k, thus sh 
abbreviates sh(k), p stands for p 1: etc. Keeping in mind that ch(y, x) = 



/ 
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ch(x, y) and sh(y,x) = sh(x, y), sh(y,x) = sh(x,y) we compute 
e B aiy x2 e- B = e B Q* XlX2 e- B = 
dyidy 2 {sh(yi, xi) ch(x 2 , y2)T> myi + ch(y h x 1 )sh(x 2 , V%)V V ^ } + 

dyidy 2 {ch(y 1 ,x 1 ) ch(x 2 ,y 2 )Q* yiy2 + sh(y u xi)sh(x 2 , y2)Q yiy2 } 

(sh o ch) (x\, x 2 ) (63a) 
3 -B _ Bq -b _ 

dy 3 dy 4 { ch(y 3 , x 1 )sh(x 2 , y 4 )V yiyz + sh(y 3 , xi)ch(x 2 , Z/^^Wj + 

dyzdy 4 { ch(y 3 ,x 1 )ch(x 2 ,y 4 )Q mm + sh(y 3 , xi)sh(x 2 , y4)Q y3y4 } 
+ (ch o sh) (xi, x 2 ) . (63b) 

Let us look first at the all the terms that come from the quartic term 
V. Since e B Ve~ 6 |0) we are interested only in terms that do not ani- 
hilate the vacuum. This implies that we keep Q* myA from (163bl) with 
everything in ( 163al) as well as Q yiy2 from ( 163al) with the last term in 
( I63bl) and the product of the two last terms in (I63ap . ( I63bl ). 

Below is the list of all the terms in j^e B Ve~ B (which do not annihilate 
|0)) in raw form, 

^ J dx 1 dx 2 dy 1 dy 2 dy 3 dy 4 ^ 

sh(yi,zi) ch(x 2 ,y2)v N {x 1 - x 2 )sh(y 3 ,x 1 )sh(x 2 ,y 4 )V myi Q* 3W + (64a) 

ch(y 1 ,x 1 )sh(x 2 ,y 2 )v N (x 1 - x 2 )sh(y 3 , x 1 )sh(x 2 ,y 4 )V yiy2 Q* my4 + (64b) 

ch(yi,xi) ch(x 2 ,y 2 )v N (xi - x 2 )sh(y 3 ,x 1 )sh(x 2 ,y 4 )Q* yiy2 Q* y3y4 + (64c) 

sh(y 1 ,x 1 )sh(x 2 ,y 2 )v N (x 1 - x 2 )sh.(y 3 ,x 1 )sh(x 2 ,y 4 )Q yiy2 Q* mm + (64d) 

(shoda)(x 1 ,x 2 )v N (x 1 - x 2 )sh(y 3 ,x 1 )sh(x 2 ,y 4 )Q* my4 + (64e) 

ch(yi,xi) ch(x 2 ,y 2 )v N (xi - x 2 ) (ch o sh) {x u x 2 )Q* yiy2 + (64f) 

(sh o ch)(xi, x 2 )v Ar (xi - x 2 )(ch o sh)(xi,a; 2 ) |. (64g) 

Next let us look at the terms that come from the cubic term, namely 
the expression e B [A, V]e _s |o). They come in two sets. First from the 
term b Xl (b xl b X2 ) there are three terms, two coming from the product of 
Q* from ( 163bl) with (162bl) and one coming from a* in f!62bj) with the 
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constant term in ( I63bl) . The are listed below 



N 1/ ' 2 / dx\dx<idy\dy<idy?,\^ 



s 



h(yi, xi)v N (xi - x 2 )(/>(x 2 )sh(?/ 2 , Xi)sh(x 2 , yz) a yi Q* y2y3 + (65a 



ch(yi,xx)v N (xi - x 2 )(p(x 2 )sh(y 2 , £i)sh(x 2 , 2/3) a* yi Q* y2V3 + (65b) 



c. 



h(yi,Xi)fAr(xi - x 2 )0(x 2 )(ch o sh)(xi,x 2 ) a* \. (65c 



The second set comes from (b* b* )b Xl and gives five terms, namely all 
terms appearing in (I63ap multiplied with a* in (162bj) . They are listed 
below, 



iV 1 I 2 i dx\dx2dyidy 2 dy3^ 



sh(yi 


xi) ch(x 2 , y 2 )v N {x x - x 2 )4>(x 2 )sh(y 3 




^ > y2yi a y 3 + 


(66a) 


ch(j/i 


Xi)sh(x 2 , y 2 )v N (x 1 - x 2 )0(x 2 )sh(y- 3 , 


Xi) 




(66b) 


ch(yi 


x-l) ch(x 2 ,y2)v N (x 1 - X2)4>(x 2 )sh(y 3 


,Xi 




(66c) 


sh(j/i 


x 1 )sh(x 2 , y 2 )v N {xi - x 2 )0(x 2 )sh(?/ 3 , 


Xi) 




(66d) 


(sh o ch) (x 1 ,x 2 )v N (x 1 - x 2 )0(x 2 )sh(y 3 , x\ 


\ * 
K 3 


}■ 


(66e) 



Some of the terms can be reduced to lower order by commuting a with 
a* whenever they appear together in a product. The term in (I66dl) 
applied to |o) gives zero. Irreducible terms are those appearing in 
fl64cl) which is fourth order, as well as the terms appearing in (I65bl) 
and f)66cp which are cubic. The quartic irreducible term in (I64cp can 
be writen, if we write ch(x, y) = S(x — y) + p(x, y) first, as follows, 

N-\l/2){ v N (y 1 -y2)sh(y 3 ,y 1 )sh(y2,y i ) + (67a) 

J dx{p(y 2 ,x)v N (y 1 - x)sh(x,y 4 )}sh(y 3 ,y 1 )+ (67b) 

J dx{p(y 1 ,x)v N (x-y 2 )sh(y 3 ,x)}sh(y 2 ,y i )+ (67c) 

dx 1 dx 2 {p(y 1 ,x 1 )p(x 2 ,y 2 )v N (x 1 - x 2 )sh(y 3 , xi)sh(x 2 , y 4 )}j . (67d) 



The rest of the quartic terms can be reduced to either quadratic or zero 
order terms. For example if we look at (I64bl) we can write, T> yim Q* y3yA = 
S(y 2 - yi) Q* yiy4 + 5{y 2 -y±) Q* yiy3 and similarly for (164a|) . Below is a list 
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of all quadatic terms, 



(1/2A0 

J 


Ucfe { 










sh(x 2 , 2/2) (sh sh) (xi, Xi)v N (xi - 


x 2 ) 


+ 


(68a) 


ch(y u x 2 ] 


sh(x ll y 2 ) (sh sh) (x h x 2 )v N (x 1 - 


x 2 ) 


+ 


(68b) 


ch(yi,ari] 


sh(x 2 , y 2 ) (sh sh) (x h x 2 )v N {x 1 - 


x 2 ) 


+ 


(68c) 


ch(yi,a;i' 


sh(xi, y 2 ) (sh sh) (x 2 , x 2 )v N (x 1 - 


x 2 ) 


+ 


(68d) 


ah(yi,xi] 


sh(x 2 , y 2 ) (sh ch) (xi,x 2 )v N (xi - 


x 2 ) 


+ 


(68e) 




ch(x 2 , y 2 ) (ch sh) (xi, x 2 )v N (xi - 


-x 2 )} . 


(68f) 



In addition (I64dl) and (|64g]) together with the trace in fl60l will supply 
a zero order term, 

\L\ = — — Tr ^(ch(A;)) 1 o m o sh(fc) + sh(fc) o m o (ch(fc)) ^ 

+ ± J dxidx2 { (sh „ ffi)^ - o *)(*,, *)+ 

(sh o sh)(xi,Xi)wAr(xi - x 2 )(sh o sh)(x 2 ,x 2 ) + 

(sh o ch)(xi,x 2 )t>Ar(xi — x 2 )(ch o sh)(xi,x 2 )| . (69) 

This term can be absorbed as a phase factor in the evolution. The 
cubic irreducible terms are those appearing in (165bj) . fl66cp . Writing 



ch(y — x) = 5(y — x) +p{y, x) we can express them in the manner below 
N~ 1/2 { v N {y 1 -y 2 )(t>{y 2 )sh{y^y 1 ) + (70a) 

'MM*-«)*«W«.*)}.h<*,») + (70b) 

dx {p(y 1 ,x)v N (x - y 2 )sh(y 3 , x)} (p(y 2 ) + (70c) 

dx x dx 2 {p(yi, x 1 )v N (x 1 - x 2 )(j)(x 2 )sh(y 2 , xi)sh(x 2 , y 3 ) } + (70e) 

dxidx 2 {p{yi, xi)p(x 2 , y 2 )v N (xi - x 2 )(fi(x 2 )sh(y 3 , xi)} } . (70f) 

The rest can be reduced to linear or eliminated, for example (166dD 
can be eliminated. For the cubic terms in the first set. From the 
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commutator relation a yi Q y2m = 5(y 1 -y 2 )a* s + 5(y 1 -y 3 )a* 2 we reduce 
them to linear and here is the list, 

N~ 1/2 J dxidx 2 { 

sh(y 1 ,x 2 ) (sh o sh) (x 1 , x 1 )4>{x 2 )v N {xi - x 2 ) + (71a) 

sh(yi, xi) (sh o sh) (xi, x 2 )(j)(x 2 )v N (xi - x 2 ) + (71b) 

ch(y 1 ,x 1 ){chosh)(x 1 ,x 2 )^(x 2 )v N (x 1 - x 2 ) \ ■ (71c) 



From the cubic terms in the second list we have the commutators 
V viy2 a l- A = <%2 - Vz)a* yi and V y2yi a* ys = 5{y x - y 3 )a* 2 hence from 
f l66al) .( |66b|) as well as from (I66el) we obtain 



N~ l/2 J dx 1 dx 2 



x 2 ) + 


(72a) 


x 2 ) + 


(72b) 


x 2 )} • 


(72c) 



6. Estimates for the error terms 
In this section we prove 
Proposition 6.1. The following estimates for the error terms holds: 

i\T 1 ||e fl Ve- B |o)|| F < CN 3 ^ 1 log 4 (l + t) 

Proof. All estimates are straightforward, based on Corollary (14. 2 p and 
Lemma (13.41) , and we only include a few typical ones from each category. 
Estimate for (I67al) : 

N^Wvsiyt - y2)sHk)(y 3 ,y 1 )sh(k)(y 2 ,y 4 )\\ L 2 (dyidy2dmdy4) 

< N- 1 \\v N \\ L o a \\ah(k)\\% < CN 30 ' 1 log 2 (l + t) 

Estimate for (167d[) : 



N 1 \\ j dx 1 dx 2 px{y 1 ,x x )p 1 {x 2 ,y 2 )v N {x 1 - x 2 )sh{y^,xi)sh{x 2 ,y A )\\ L ^ dyidy2dy3dy ^ 

< iV- 1 ||^|| L oo||sh(A;)||i 2 ||p 1 ||i 2 < CN 3 ^ 1 log 4 (l + t) 
Estimate for f!68fj) . keeping only the 5 contribution in ch(fc): 
N- 1 \\sh{k){ yi , y 2 )v N { yi ~y 2 )\\ L 2 {dyidy2) < \\v N \\ L °o\\sh(k)\\ L i < CN 3 ^ 1 log(l + 1) 
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Estimate for (170 aft . : 

N~ 1,2 \\v N (y 1 - y 2 )(f){y2)sh{k){y^yi)\\ L 2 {dyid y 2dy3) 

< |MMMMM*)I|l» < CN^ hgi ^ f) 

Estimates for (172bj) : 



N 1/2 || / dx 1 dx 2 p 2 (y 2 ,X 2 )(sh(k) Osh(k))(x 1 ,X 1 )(j)(x 2 )v N (x 1 - X 2 )\\ L 2^ dy2 



< N 1/2 ||p 2 ||L 2 |l sri (^) ° sh(A;)(a;,a;)||i,i( d:r )||0||i,oo||u J v||L2 
<iV- 1 / 2 ||j9 2 || L2 ||sh(A;)|| 2 :2 ||0|| L o O ||t; Ar || L2 

<Civ^ log3(1 + t) 

- OJV f 3/2 

□ 

7. Appendix 

The purpose of this appendix is to make certain connections with 
our previous work [20], [21] and to provide some useful information. 
We would like to explain first the Lie algebra isomorphism which is 
crucial in our work. This was explained in |20] but we include it here 
for completeness and for the convenience of the reader. Let us define 
first 



and using these form 

A(f) := J dx{A T J{x)} = J dx{Mx)a x + f 2 (x)a* x } . (73) 
It is straightforward to check the commutation, 

[.4(f), .4(g)] = J dx{f 1 (x)g 2 (x)-f 2 (x)g 1 (x)} 

= - J dxdy {f T (x)S(x - y) My)} . (74) 

For L 2 (dxdy) kernel functions d(t,x,y) and k(t,x,y), l(t,x,y) such 
that k and I are symmetric in the (x, y) variables form the symplectic 
matrix 

S(d,kJ):= ( d _ k ^) . (75) 
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We will write S(x,y) when convenient, suppressing the t dependence. 
Next we define the map X : sp(C) (-> Quad from the space of complex, 
L 2 symplectic matrices to quadratic expressions in (a, a*) as follows : 

X{S) := ~ J dxdy{AlS(x,y)JA y } . (76) 

Theorem 7.1. Let f(x) a vector function and A(f) the corresponding 
expression (173]) . Then the following commutations relations hold 

[l(S),A(f)]=A(Sof) (77a) 

e I ( s U(f)e- I ( s ) = ^(e s of ) (77b) 

where (for example) Sot = J dx{S(x,y)f(y)} etc. Formula ( 177b|) holds 
for any complex symplectic S G sp c (lR) = (C T ) _1 sp(lR)C T where C is 
the change- of -basis matrix 

1 (I -il 

and sp(lR) is the Lie algebra of real symplectic matrices. This condition 
insures that e s is unitary. 

Proof. The commutation relation in (I77ap can be easily checked, but 
we point that it follows from (!7ij) . For any rank one quadratics we 
have using (!7i|) 

[A(f)A(g) , A(h)} = -A((gf T + fg T ) o Jh) . 
Thus for any R we have 

J dxdy{A T x R{x,y)Ay} , A(f) = -A((R + R t ) o Jf ) . 

Now specialize to R = {1/2)SJ with S G sp and use S T = JSJ to 
complete the proof. 

For the second formula, introduce a complex parameter t and take 
\ip) a Fock space vector with finitely many non-zero components. It is 
trivial to check, using flTTap, that all t derivatives of e tx ^ A{f)e- tx ^\i)) 
and A{e ts of ) |^) agree when t = 0, and both the left hand side and 
right hand side are analytic if \t\ is sufficiently small. Thus they agree 
for all t complex, sufficiently small. To take t large, we have to restrict 
ourselves to real t and use the group properties of the unitary family 



e ts . 



A more formal (but convincing) "proof follows from 

1 
2! 



e T(s) Ce -x(s) =c+ [X{S)A + L\Z(S), [X{S),C\] + 
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□ 

Theorem 7.2. The map X : sp(C) (->■ Quad defined in (175]) a Lie 

algebra isomorphism. Moreover for S, R G sp c (R) we have the formulas 

l[(d t e s ) o e - 5 ) = (ft^e- 1 ^ (78a) 

j(e 5 o R o e' 5 ) = e x ^l(R)e~ x ^ . (78b) 

Proof. We point out that (177bl) implies (I78bl) . at least when R is rank 
one, R := f(x)g T (y). Notice that (177bl) can be written 

e 1 ^ J dx{f T (x)A x }e~ x ^ = J dxdy [f T (x)e sT Ay) 
from which we have 

J dxdy{A T J{x)% T {y)JA y }e- x ^ 

= <FWA(f)e- x W f dy{g T (y)JA y }e- x ^ 

= A{e s f) J dydz{g T (y)Je JSJ A z } 

= J dxdy{A T 'e s Re~ s JAy} 

since S T = JSJ and Je JSJ = e~ s J. 

A direct proof of the Lie algebra isomorphism follows from the fol- 
lowing elementary computations. Define first 

Qxy ■= a x a y and Q* xy := a* x a* and X xy := -{a x a* + a* y a x ) . 
It is straightforward to verify, 

[Q X1X2 , Q* yiy2 ] =S(x x - yi)M X2y2 + 5(x 1 - y 2 )M X2yi + 5(x 2 - yi)M xiy2 + 5(x 2 - y 2 )M x 
[Qx lX2 ,N yiy2 ] =8{x x - y 2 )Q X2Vl + 5(x 2 - y 2 )Q xiyi 

[4, 2 > Q* yiy2 ] = - yi)Q* X2y2 + Km - yi)Ql 2yi 

[N XlX2 ,N yiy2 ] =5(xi- y 2 )M yiX2 -5(x 2 - yi)M XlV2 
and using these we can directly verify that 

i J dx 1 dx 2 {k(x 1 ,x 2 )Q XlX2 } , - J dy 1 dy 2 {l(y 1 ,y 2 )Q* yiy2 } 
= - J dxdy{(kol)(x,y)M xy } 
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which corresponds to the relation 



k 






1 



kol 
-I ok 



The other cases are checked in a similar manner. 

To prove f)78al) . expand both the left and right-hand side as 

l((d t e s )e- s )=l(s + 1 - [ S,S] + ..>) 
= ±(S) + ~[l(S),i(S)] +... = {d t e^ s) )e-^ s) . 



The proof of (178bj) is along the same lines. The proofs can be made 
rigorous by an analyticity argument on the dense subset of vectors with 
finitely many non-zero components. □ 



The second part of this appendix is devoted to the connection be- 
tween the present equations for the pair excitations ( 117bl) . (I17cf) with 
the evolution equation derived in [20] and [21], the idea being that the 
nonlinear equation derived in these references turns out to be linear 
when expressed in new coordinates. The connection however is not 
entirely obvious. From the identities, 



2K = (ch{2k) sh(2£Q \ 
" v sh(2A;) ch(2A;) J 

( ch(fc) o ch(A;) + shpO o sh(A;) ch(A;) o sh(/c) + sh(/c) o ch(k)\ 
' \sh(k) o ch(A;) + ch(fc) o sh(A;) sh(A;) o sh(fc) + ch(fc) o ch(fc) / 

_ /2sh(fc) osh(fc) + 1 2sh(A;) o ch(fc) \ 
V 2ch(fc) o sh(fc) 2sh(fc) o sh(fc) + l) 

and 



/ ch(fc) o ch(fc) - sh(k) o sh(k) -ch(k) o sh(k) + sh(k) o ch(A;)\ 
\sh(k) o ch(fc) - ch(A;) o sh(k) -sh(k) o sh(k) + ch(k) o ch(A;) / 



we have ch(2A;) = 2sh(A;) o sh(k) + 1 = 2ch(fc) o ch(k) - 1 and sh(2fc) = 
2ch(A;) o sh(k) = 2sh(fc) o ch(A;). Using the above identities we can 
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readily derive the identities, 



W (sh(Jfe) o sh(A;)J = S (sh(Jfe)) o sh(k) - sh(Jfe) o S (sh(Jfe)) (79) 
= W fch(fc) o c h(k)) = W (ck(k)) o ch(fc) + ch(fc) o W (ch(k) 



S (ch(k) o sh(ife)) = W (ch(fc)) o sh(Jfe) + ch(fc) o S (sh(fc)) (80) 

= S (sh(Jfe) o ch(A;)) = S (sh(Ar)) o ch(k) - sh(jfe) o W (ch(k) 

The equations above have a generic character when we consider for 
example S(/ o g) = W(J) o g + J o S(#) etc. 
Our original evolution equation, see [21] was 

-sh(k) t + g T o sh(fc) + sh(fc) o <? — ch(&) o m^j o ch(/c) = (81) 

-ch(fc) 4 + [g T , ch(A;)] + sh(&) o o sh(A;) , 
or its equivalent form 

1.,.. T 



sh(fc)t + </ J o sh(fc) + sh(&) oj- ch(k) o m (82) 

^ch(fc) 4 + [# T , ch(fc)] + sh(fc) o mj o (ch(fc)) -1 o sh(Jfe) = . 
The second equation can be written, 

S(sh(fc)) - W(ch(fc)) o (chikjy 1 osh(k) = ch(fc) o m + sh(fc) o m o (ch(fc)) _1 sh(fc) . 
Multiplying on the left ([82]) with (ch(A;)) 1 we obtain 

(ch(fc)) _1 oS(sh(fc)) - (ch(A0) _1 o W(ch(T)) o (chtT)) -1 o sh(fc) = 

m + (ch(fc)) sh(fc) o mo (ch(Jfe)) osh(A;) 
and using the fact that 

(clip0) _1 W(^clip0)(cR^)) _1 = -w((ch(fc)) _1 ) 
we have in view of the identities ( 180]) the simple equation, 

S(() = m + (omo(, (83) 

where ( := (ch(/c)) 1 o sh(k) = sh(k) o fch(A;)) . Now from (1831 we 
can readily derive the equation below, 



W(sh(/c) o sh(/c)) = mo sh(fc) o ch(fc) — sh(fc) o ch(fc) o m (84) 
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which is the same as ( I43bl) and it already appeared in [2T]. To derive 
an equation for sh(2fc) notice that we can write, 

sh(2A;) = (ch(F)) 2 C + C(ch(A;)) 2 

and subsequently compute using fl80|) .f l83|) and fl84|) in a tedious but 
straightforward manner, 

S(sh(2*)) = W((ch(fc)) 2 )C + (du^)) 2 S(C) + S(C)(ch(A;)) 2 - CW((ch(fc)) 2 ) 
= \m o sh(/c) o ch(/c) — sh(/c) o ch(/c) o m] o (ch(fc)) -1 o sh(/c) 
+ (ch(A;)) 2 o m + (ch(fc)) 2 o (ch(fc))- 1 o sh(fc) o to o (ch(fc)) _1 sh(fc) 
+ m o (ch(A;)) 2 + sh(fc) o (ch(A;))~ 1 o m o sh(fc) o (ch(A;))~ 1 o (ch(A;)) 2 
— sh(fc) o (ch(fc)) o [to o sh(fc) o ch(fc) — ch(fc) o sh(fc) o ml 
= m o ch(2A;) + ch(2A;) o m . 

This equation is ( 143 aft . 
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